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Abstract 

The design of strategies to generate efficient mixing is crucial for a variety of ap- 
plications, particularly digital microfluidic devices that use small "discrete" fluid 
volumes (droplets) as fluid carriers and microreactors. In recent work, we have 
presented an approach for the generation and control of mixing inside a translat- 
ing spherical droplet. This was accomplished by considering Stokes' flow within a 
droplet proceeding downstream to which we have superimposed time dependent (si- 
nusoidal) rotation. The mixing obtained is the result of the stretching and folding 
of material lines which increase exponentially the surface contact between reagents. 
The mixing strategy relies on the generation of resonances between the steady and 
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the unsteady part of the flow, which is achieved by tuning the parameters of the 
periodic rotation. Such resonances, in our system, offer the possibihty of controUing 
both the location and the size of the mixing region within the droplet, which may 
be useful to manufacture inhomogeneous particles (such as Janus particles). While 
the period and amplitude of the periodic rotation play a major role, it is shown 
here by using a triangular function that the particular shape of the rotation (as a 
function of time) has a minor influence. This finding demonstrates the robustness 
of the proposed mixing strategy, a crucial point for its experimental realization. 

Key words: Digital Microfluidics, Droplet, Stokes' flow. Chaotic mixing. 
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1 Introduction 

Although most microfluidics have been using fluid-streams as the main means 
to carry fluids, devices based on individual droplets have been proposed as 
well. In the latter, also called "digital microfluidic" systems, "discrete" fluid 
volumes (droplets) rather than continuous streams are used, with the poten- 
tial to utilize individual droplets as microreactors [1]. 

Whether fluids are encapsulated within droplets or flow along channels, re- 
actions can occur efficiently only in presence of rapid mixing. Such mixing 
conditions are not easy to fulfill due to the low Reynolds number of the fiow 
involved, which prevents turbulence from taking place. Stirring, in addition 
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to molecular diffusion, is needed to stretch and fold fluid elements, thus sig- 
nificantly increasing interfacial areas. Whereas some strategies are based on 
complex channel geometry for flows in microchannels, active methods (via 
external forcing) (see, e.g. [2ll3|l¥|l5|l6] ) have also resulted in efficient mixing, es- 
pecially at very low Reynolds numbers [7j. The combination of both geometry 
alteration and forcing has been explored as well [TfSfQllTOlITT] . 
Chaotic advection inside a liquid droplet subjected to a forcing has been stud- 
ied extensively |12|ll3fl4fl5|ll6fl7fl8j and obtained experimentally by means of 
oscillatory flows p^pO] . In this paper as in pi], we concentrate on controlling 
both the size and the location of the mixing. Indeed, we are interested in both 
complete and incomplete mixing within drops. While complete mixing should 
be useful for efficient reactions to occur uniformly, there are applications in 
which non-uniform mixing, when well-controlled, could be desired. The latter 
include the formation of inhomogeneous drops and particles, such as Janus 
particles, whose properties have been shown, in certain cases, to be superior 
to their homogeneous counterparts. In this paper, we extend the results of [21] 
which were obtained by using unsteady - sinusoidal - forcing on a translating 
droplet. Specifically, we study the robustness of the control of the size and 
location of the mixing region by selecting a different time function for the 
rotation. 

The existence of chaotic behavior in three-dimensional bounded steady flows 
has been shown (e.g. p^l3f23|l2l] ). Our work is distinct from the latter con- 
tributions through the addition of unsteadiness, which is crucial to control the 
chaotic mixing behavior through resonance effects [25]|26|I18] . 
The physical system under consideration in this work is described in Sec. |2} 
and our numerical results are displayed in Sec. |3j The control of both the 



location (Sec. 3.1) and size (Sec. 3.2) of the chaotic mixing region is studied 
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qualitatively via Liouvillian sections. 



2 Physical system 

2.1 Flow equations and assumptions 

Let us assume a Newtonian liquid droplet of sperical shape suspended in an 
incompressible Newtonian fluid. Its motion consists of a translation and a slow 
rigid body rotation (see y^). As in the previous reference, we assume that the 
interfacial tension is sufficiently large for the drop to remain spherical through 
its motion and that the Reynolds number is very small compared to 1. Thus, 
a reasonable approximation is to consider that both the internal and external 
flows are Stokes flows. The boundary conditions at the droplet surface can be 
derived from the continuity of velocity and tangential stress conditions. 
The resulting internal flow is a superimposition of a steady base flow (non- 
mixing flow) and an unsteady rigid-body rotation. Consider a Cartesian coor- 
dinate system translating with the center-of-mass velocity of the droplet with 
the orientation of the axes such that the unit vector points in the direction 
of the translation and the unit vector lies in the a; — plane. We then 
obtain the following unsteady internal flow: 

u = zx- ea^{t)u^y, (1) 
v = zy + ea^it) {u];,x - uJ3:z) , (2) 




where v and w are the time derivatives of x, y and z, respectively. In 
Eqs. ([l]|3|, all the lengths and velocities were made dimensionless by using the 
droplet radius and the magnitude of the translational velocity as the length 




(3) 
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and velocity scales, respectively. The rotation is characterized by a rotation 
with a maximum amplitude e -C 1, a fixed orientation vector 

uj = (v^/2, 0, , 

and a periodic amplitude of the rotation a^(t). In the present paper we study 
two time dependent functions for the amplitude of the rotation a^^lt), one 
consisting of only one harmonic: 

= ^(l + cosu;t), (4) 

and the other one consisting of a triangular function with an infinite number 
of harmonics: 

%{t) = 1 + ^a2n^icos{{2n - l)ujt) , where a„ = f — j . (5) 
n=i Kmc J 

Equations ([l](3]) are the same as in [13] except that the vorticity vector, that 
was constant in [13], is unsteady in our case. Note that the internal fiow is a 
solution of the (time dependent) Stokes problem in which a time-dependent 
body force has been added on the right hand side of the momentum equation. 
In practice, this time-dependent forcing could be realized, e.g. by creating a 
time-dependent swirl motion in the external fiow or by applying an electric 
field that exerts a torque on the drop (see, e.g. [27] or work on electrorotation 
[28]). Note that the surface of the droplet, = + + 2;^ = 1 is invariant 
under fiow (|T]|3|). 

2.2 Non-mixing case 

The non-mixing (unperturbed) case or base fiow, i.e. e = 0, is characterized 
by two invariants (i.e. time-independent quantity), the stream function ip and 
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the azimuthal angle such that: 



= 2 (^^ ~^ ~ ^'^ ~ y"^ ~ -^^) ' ^ ~ arctan {y/x) , 

where ip E [0, 1/8] and (p ^ [0,27r[. The streamhnes of the non-mixing case are 
joint hnes of constant ip and 0, denoted by T^^^. The heterochnic orbits T^^^ 
such that ip = and G [0,27r[, connect two hyperbohc fixed points located 
at the poles of the sphere (see Fig. [T]). All other streamlines are closed curves 
that converge toward a circle of degenerate elliptic fixed point (x^ + = 
1/2,2; = 0) as '0 is increased toward the value 1/8. The frequency of motion 
on the streamline ^ is independent of and given by 

i-v2 y^i + 7(^) sin/3 vVi + 7y' 

with 7('0) = Vl — 8'0 and K is the complete elliptic function of the first 
kind. The plot of as a function of t/', presented in Fig. [ij shows that is 
bounded by two limits, i7(0) = and i7(l/8) = On every streamline F^^,^ 
we introduce a uniform phase x mod 27r such that x = on the — e.y plane 
and X = ^ (V")- With this uniform phase, the non- mixing flow can be described 
by using the variables (-0, 0, x) instead of the Cartesian coordinates (x, z): 

^- = 0, = 0, x = ^^(V')- 

Such a system is generally classified as action- action- angle ^ where the two 
invariants -0 and are the two actions and x is the angle. 
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Fig. 1. Streamlines in a cross-section of the drop (without rotation) and their 



2.3 Mixing case 

In the mixing (or weakly perturbed) case < e ^ 1, the quantities ip and (p 
are no longer constant and their time evolution satisfies 



where G{ip,x) = z/{x'^ + y^) is 2tt periodic in x ^^id has zero average in x- 
The time evolution equation for x reads 



where H is 27r periodic in x- Such a system possesses two time scales, a fast 
one (of order one) associated with x, typically on time scales of 1, and a slow 
one (of order 1/e) associated with both ip and 0. The generation of chaos 
arises from the resonances between the frequency of the integrable case Vt{il)) 
and the forcing frequency uj. 



frequencies {ip) as given by Eq. 



ip = -2ea^{t)u^ij sin 0G {ip, x) , 



(p = eauj{t)u^ - ea^{t)u^ cos (pG {ip, x) , 



X = n{ij) + ea^{t)H{ij,^,x), 
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3 Generation and control of the chaotic mixing region 

In this work we fallow the approach of [21] that enables the generation of a 
three-dimensional chaotic mixing region, for which we are able to control both 
the location and the size. The method consists in bringing a specific family 
of the unperturbed tori { \ into resonance with the periodic per- 

turbation a^{t) by selecting the frequency u so that it satisfies the resonance 
condition: 

nfl {4'n,m) — = 0, for n, m G N^. 

The control is realized by selecting the two parameters that characterize the 
periodic rigid body rotation, specifically its the maximum amplitude e and its 
frequency uj. The amplitude satisfies < e << 1, where the lower and upper 
limits correspond to the absence of mixing and complete mixing, respectively. 
In this work we are interested in the location of the mixing regions within 
the drop. To study our system, we compute two-dimensional projections of 
time-periodic three-dimensional fiows via the combination of a stroboscopic 
map and a plane section (here, the y = plane). In other words, the points 
on these so called Liouivillian sections are the intersections of the trajectories 
with the plane ?/ = at every period 2ti/uj. 

3. 1 Control of the location 

In this paper, we study the infiuence on mixing for two different time- 
dependent functions of the periodic forcing: C^^it) that corresponds to only 
one harmonic and T^{t) that is a triangular shaped forcing (see Fig. [2] and 
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Fig. 2. Right panel: Spherical droplet represented with rigid body rotation (axis 
in dash line). Middle panel: Plot of the periodic forcing T^{t) of the rigid body 
rotation. Left panel: Plot of the periodic forcing Cuj{t) of the rigid body rotation. 



Figure [3] shows the case where the periodic forcing a^^it) of the rigid body ro- 
tation has only one harmonic at^(t) = Here, one can observe two major 
three-dimensional chaotic mixing regions, all others being of almost negligible 
size. The first one of these major mixing regions is around the pole-to-pole 
connections in the center and near the surface of the droplet. This region 
corresponds to the resonance of the family of unperturbed tori |r^^ ^| with 
n > 1 and is labeled as the central chaotic mixing region. The other one 
corresponds to the resonance of the unperturbed tori F^^^ and is labeled as 
the primary chaotic mixing region. While the central mixing region is always 
present at this particular location inside the drop, the primary one can be 
displaced within the droplet. For small values of cj, the primary and central 
chaotic mixing regions coincide, located around the pole-to-pole connection. 
For larger values of cu, the primary chaotic mixing region separates from the 
central one, and penetrates deeper into the droplet. As uo is increased further, 
it moves toward the location of the circle of degenerate elliptic fixed point of 



Eqs. dip] )). The two functions differ by: 



T^W-CWI <0.11, VtGM. 
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the unperturbed system {z = 0, + = 1/2) by following the location 
given by ipi^i. 

Figure |4] shows the case where the periodic forcing of the rigid body rotation 
is ai^(t) = %o{t). As in the situation where only one harmonic is present, one 
observes two major chaotic mixing regions. While the central region remains 
around the pole-to-pole connection, the primary one can be shifted across the 
droplet by varying the frequency u, as it is the case for C^(t). 
Although the triangular amplitude 7^ is composed by an infinite number of 
harmonics, only the first harmonic plays a non-negligible role. All other chaot- 
ic mixing zones created by higher harmonic appear of negligible size within 
the values of parameters studied. This result reinforces the robustness of the 
proposed method and confirms that the main parameter influencing the loca- 
tion of the primary mixing zone is the frequency u, and not the specific shape 
of the forcing. 



3.2 Control of the size 



Whereas the tuning of the frequency uj of the rigid body rotation enables the 
placement of the primary chaotic mixing region at a particular location inside 
the droplet, it is its amplitude e that mostly determines the size of the main 
mixing regions. 

The sizes of the primary and central chaotic mixing regions are presented in 
Figs. 5][6 for the two perturbations discussed above. It is clear that in both 
cases the size increases with the amplitude of the rotation e. Such a control 
of the size of the mixing allows one to vary the extent of mixing generated 
within the droplet from no mixing to complete mixing. 
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Fig. 3. (Color on line) Liouvillian sections for the time-dependence C^^ with the 
frequencies uj = 0.95,1.10,1.25,1.40 (a-d) and the amplitude e = 0.05. The (red) 
dashed line inside the primary mixing region is T^^ ^ . 

Once again it is important to note that for the case where the periodic amph- 
tude of rotation is triangular only the first harmonic plays a significant role. 
Figure |6] shows the case of nearly no mixing ( e = 0.01) and that of complete 
mixing (e = 0.20), the predominance of the first harmonic is not affected by 
the magnitude of e. 

4 Conclusions 

In this paper, we have shown that a chaotic mixing region within a droplet is 
generated by adding a slow oscillatory rigid-body rotation to the translating 
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Fig. 4. (color online) Liouvillian sections for the time-dependence 7^ with the 
frequencies oo = 0.95, 1.10, 1.25, 1.40 (a-d) and the amplitude e = 0.05 with 
o-wit) = %jit)- The (red) dashed line inside the primary mixing region is the torus 

drop motion. Our strategy consists in using resonances between one of the 
frequencies of the base flow generated by the steady translation and the forcing 
frequency characteristic of the oscillatory rigid body rotation. This strategy 
allows for a direct control of both the location and the size of the mixing by 
judiciously adjusting the parameters of the time- dependent rotation, i.e. its 
amplitude and frequency. One the one hand, controlling the mixing location 
could be useful to manufacture inhomogeneous droplets and particles (such as 
Janus particles) whose properties could be superior to homogeneous ones. On 
the other hand, controlling the size of the mixing zone could be important to 
control reaction rates. Note that in this system, the total mixing state does 
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Fig. 5. (Color on line) Liouvillian sections for the time-dependence C^^ with the 
frequencies uj = 1.25 and the amplitude e = 0.01,0.05,0.10,0.20 (a-d). The (red) 
dashed line inside the primary mixing region is T^^ ^ . 

not have islands of non-mixing zones, as is often the case with other mixing 
strategies. Size and location control was obtained not only with a sinusoidal 
rotation (with only one harmonic) but also with a rotation consisting of a time 
dependent triangular function, with an infinite number of harmonics. In the 
latter case, the first harmonic only was shown to play a significant role. 

Acknowledgements 

This article is based upon work partially supported by the NSF (grants CTS- 
0626070 (N.A.), CTS-0626123 (P.S.) and 0400370 (D.V.)). D.V. is grateful 



13 



Fig. 6. (Color on line) Liouvillian sections for the time-dependence with the 
frequencies uj = 1.25 and the amplitude e = 0.01,0.05,0.10,0.20 (a-d). The (red) 
dashed line inside the primary mixing region is T^^ ^ . 

to the RBRF (grant 06-01-00117) and to the Donors of the ACS Petroleum 
Research Fund. CC acknowledges support from Euratom-CEA (contract 
EUR 344-88-1 FUA F) and CNRS (PICS program). 

References 

[1] H. Song, J.D. Tice and R. F. Ismagilov, A Microfluidic System for Controlling 
Reaction Networks in Time, Angew. Chem. Int. Edit. 42, (2003) 768-772. 

[2] M.H. Oddy, J.G. Santiago, J.C. Mikkelsen, Electrokinetic Instabihty 
Micromixing, Anal. Chem. 73, (2001) 5822-5832. 



14 



[3] H.H. Bau, J. Zhong and M. Yi, A minute magneto hydrodynamic (MHD) mixer, 
Sens. Actuators B 79, (2001) 207-215. 

[4] A. Ould El Moctar, N. Aubry and J. Batton, Electro-hydrodynamic micro- 
fiuidic mixer. Lab Chip 3, (2003) 273-280. 

[5] I.K. Glasgow and N. Aubry, Enhancement of microfiuidic mixing using time 
pulsing. Lab Chip 3, (2003) 114-120. 

[6] I.K. Glasgow, J. Batton and N. Aubry, Electroosmotic mixing in microchannels. 
Lab Chip 4, (2004) 558-562. 

[7] A. GouUet, I.K. Glasgow and N. Aubry, Effects of microchannel geometry on 
pulsed flow mixing, Mech. Res. Commun. 33, (2006) 739-746. 

[8] X. Niu and Y-K. Lee, Efficient spatial-temporal chaotic mixing in 
microchannels, J. Micromech. Microeng. 13, (2003) 454^462. 

[9] F. Bottausci et ai, Mixing in the shear superposition micromixer: three- 
dimensional analysis, Phil. Trans. Royal Soc. A 362, (2004) 1001-1018. 

[10] M.A. Stremler, F.R. Haselton and H. Aref, Designing for chaos: applications 
of chaotic advection at the microscale, Phil. Trans. Royal Soc. A 362, (2004) 
1019-1036. 

[11] Y.K. Lee, C. Shih, P. Tabeling, CM. Ho, Experimental study and nonlinear 
dynamic analysis of time-perioidc micro chaotic mixers, J. Fluid Mech. 575, 
(2007) 425-448. 

[12] K. Bajer and H.K. Moffatt, On a class of steady confined Stokes flows with 
chaotic streamlines, J. Fluid Mech. 212, (1990) 337-363. 

[13] D. Kroujiline and H.A. Stone, Chaotic streamlines in steady bounded three- 
dimensional Stokes flows, Physica D 130, (1999) 105-132. 



15 



[14] S.M. Lee, D.J. Im and I.S. Kang, Circulating flows inside a drop under time- 
periodic nonuniform electric fields, Phys. Fluids 12, (2000) 1899-1910. 

[15] T. Ward and G.M. Homsy, Electrohydrodynamically driven chaotic mixing in 
a. translating drop, Phys. Fluids 13, (2001) 3521-3525. 

[16] R.O. Grigoriev, Chaotic mixing in thermocapillary-driven microdroplets, Phys. 
Fluids 17, (2005) 033601.1-033601.8. 

[17] X.M. Xu and G.M. Homsy, Three-dimensional chaotic mixing inside drops 
driven by a transient electric field, Phys. Fluids 19, (2007) 013102.1-013102.11. 

[18] D. Vainchtein, J. Widloski and R. Grigoriev, Resonant chaotic mixing in a 
cellular flow, Phys. Rev. Lett. 99, (2007) 094501.1-094501.4. 

[19] T. Ward and G.M. Homsy, Electrohydrodynamically driven chaotic mixing in 
a translating drop part II: Experiments, Phys. Fluids 15, (2003) 2987-2994. 

[20] R.O. Grigoriev, M.F. Schatz and V. Sharma, Optically controlled mixing in 
microdroplets. Lab Chip 6, (2006) 1369-1372. 

[21] R. Chabreyrie et al, Tailored mixing inside a translating droplet, Phys. Rev. E 
77, (2008) 036314.1-036314.4. 

[22] M. Feingold, L.P. Kadanoff and O. Piro, Passive scalars, three-dimensional 
volume-preserving maps, and chaos, J. Stat. Phys. 50, (1988) 529-565. 

[23] D. Vainchtein, A. Vasiliev and A. Neishtadt, Adiabatic chaos in a two- 
dimensional mapping. Chaos 6, (1996) 514-518. 

[24] D. Vainchtein, A. Neishtadt and I. Mezic, On passage through resonances in 
volume-preserving systems. Chaos 16, (2006) 043123.1-043123.11. 

[25] R. Lima and M. Pettini, Suppression of chaos by resonant parametric 
perturbations, Phys. Rev. A 41, (1990) 726-733. 



16 



[26] J.H.E. Cartwright, M. Feingold and O. Piro, Chaotic advection in three- 
dimensional unsteady incompressible laminar flow, J. Fluid Mech. 316, (1996) 
259-284. 

[27] N. Aubry and P. Singh, Influence of particle-particle interactions and particle 
rotational motions in traveling wave dielectrophoresis. Electrophoresis 27, 
(2006) 703-715. 

[28] W.M. Arnold, U. Zimmermann, Electro-rotation: development of a technique 
for dielectric measurements on individual cells and particles, J. electrostat. 21, 
(1988) 151-191. 



17 



